function plot_r()
time = dlmread('depos_MPC_t');
r2 = dlmread('depos_MPC_r2');
plot(time,r2);
set(gca,'FontSize',24);
xlabel('time(s)');
ylabel('<r^2>');
